Excitations in one-dimensional S— 1/2 quantum antiferromagnets 
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The transition from dimerized to uniform phases is studied in terms of spectral weights for spin 
chains using continuous unitary transformations (CUTs). The spectral weights in the 5 = 1 channel 
are computed perturbatively around the limit of strong dimerization. We find that the spectral 
weight is concentrated mainly in the subspaces with a small number of elementary triplets (triplons), 
even for vanishing dimerization. So, besides spinons, triplons may be used as elementary excitations 
in spin chains. We conclude that there is no necessity to use fractional excitations in low-dimensional, 
undoped or doped quantum antiferromagnets. 
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] The determination of elementary excitations is a fun- 
damental issue in condensed matter physics. Especially 
in the context of high-temperature superconductivity, the 
role of magnetic excitations in two-dimensional quantum 
antiferromagnets is strongly debated 0, S S 01 ■ The 
bottom line of the debate is whether or not it is possi- 
ble and/or necessary to consider fractional excitations, 
so-called spinons with S = 1/2, or whether integer exci- 
tations, magnons or elementary triplets, triplons [2J, with 
5 = 1, can equally be used. Here we provide clear indica- 
tion that even in one-dimensional (ID) spin chains where 
spinons are a valid description an alternative description 
by triplons is possible. Hence the necessity to resort to 
fractional excitations is quite generally questioned. 

The spinon concept comes from ID spin chains (5 = 0) 



H' = J £(1 + *(-!)*) SiS 4 



(1) 



where spinons with S = 1/2 are established as the el- 
ementary excitations 0, la]- For any finite dimeriza- 
tion 5^0, however, the spin chain can be described 
by massive triplets with S = 1. They can be viewed 
as bound states of two spinons 0, 0, 0, 0] which 
are, however, not free due to their confining interac- 
tion. Starting from strong dimerization, where a de- 
scription with triplons is certainly appropriate, the ques- 
tion arises what happens to the triplons for vanishing 
dimerization. Besides its academic interest, this ques- 
tion is of great practical importance since in recent years 
triplon-based techniques have been advanced consider- 
ably H E E E 0. Our recent progress 
has rendered even the computation of spectral densities 
possible 

In this work, we investigate to which extent a triplon- 
based description is possible by computing spectral 
weights of multi-triplon contributions for dimerized spin 
chains. Two scenarios are conceivable for S — > 0: 
(i) The spectral weight is distributed rather evenly over 
the multi-triplon channels implying that the weight of a 



particular channel is small and that a large number of 
channels must be taken into account, see discussion in 
Ref.El This scenario would make an approach in terms 
of triplons difficult to use, hence inappropriate, 
(ii) The spectral weight is found mainly in the channels 
with a small number of triplons, implying that a triplon 
approach is very useful and appropriate because spectral 
properties can be computed from the dynamics of a small 
number of excitations in a small number of channels. 

Scenario (i) Let us look at a field theoretic descrip- 
tion for indications for this scenario. Using abelian 
bosonization 0, Eii. E3. l2fij the low-lying states of an 
anisotropic spin chain are described by the Hamiltonian 



Hft — 



2- 



[K(wU(x)) 2 + K- l {d x ${x)Y] alx (2) 



where v cx Jo is the spin-wave velocity and K is the 
interaction parameter, e.g. K = 1 for the XY-model, 
equivalent to free fermions, and K = 1/2 for the isotropic 
case. A local operator is 

S* = (2n)- 1 d x 'i>(j)+A(-iyco S (2^(j)) (3) 

where the lattice constant is unity, and A a non-universal 
constant. The undimerized Hamiltonian Q corresponds 
essentially to free bosons (J2J. A finite dimerization 
is accounted for by an additional term proportional to 
5 cos(2$(x))g?x leading to a sine-Gordon model. For 
K < 2 the system is massive with a gap A cx (jVt 2 -^) 
(= <5 2 / 3 for the isotropic chain). 
A single bosonic mode of is created by 



61 = 1 



2 $(k)/N k - tN k IL(k) 



(4) 



where N k = \JixKj\k\. The quantities with tilde are the 
Fourier transforms of the real-space fields. We expand 
the excited state SJ|0) in states of various number of 
bosons with focus on the vicinity of momentum 7r where 
most of the weight in the dynamic structure factor is 



2 



found: oc cos(2$(j)). The coefficients of one- and 
two-boson states are 

c fc = <0|6 fc cos(2$)|0) = (0|[6 fc ,coB(2$)]|0) (5a) 
c k , q = (0|6 fc 6 g cos(2$)|0) = <0|[fo fc [6 g ,cos(2$)]]|0)(5b) 



We make use of 



[b k ,2$(j)] = s/K/W\V^/Lexp(ikj) 



(6) 



where the momenta are discrctized by a finite system size 
L to ensure normalizability. This yields c k = and 



Ck,q OC 



K 2tt 



exp(i(k + q)j)A K 



(J) 



where (sin(2<f>)) = and (cos(2$)) cx A K is used. So the 
total weight W* in the one-boson channel is always zero. 
The total weight W£ in the two-boson channel is 



oc A 2K ln(A) 2 
oc S 2K ^ 2 - K Hn(S) 2 



(8) 



in leading order in ln(<5). Generally, all channels with an 
odd number of bosons carry no weight (at momentum it) 
whereas channels with 2n bosons have 



WL oc 5 2K '^ 



ln(<5) 5 



(9) 



whence we conclude that any single channel becomes ir- 
relevant on 5 — ► 0. Only the consideration of an infinite 
number allows to treat vanishing dimerization correctly. 
This appears to be sound evidence for scenario (i). 

Scenario (ii) Our main argument is an explicit cal- 
culation of spectral weights for the first triplon channels. 
First, however, we challenge the general validity of the 
previous argument by four basic considerations. 

(a) Considering free fermions (K — 1) we know that 



the operator Sf 



1/2 excites a particle-hole contin- 



uum at all wave-vectors. So the dynamic correlations of 
this operator are exhaustively described by two elemen- 
tary excitations, a particle and a hole, independent from 
the dimerization 5. 

(b) Considering the isotropic spin chain (K — 1/2) it 
was shown that 72.89 % of the total weight (sum over all 
wave vectors) of the dynamic structure factor is described 
by the fwo-spinon continuum |26| . 

(c) It is undisputed that the dimerized isotropic spin 
chain displays a single-mode peak with finite spectral 
weight at all wave vectors as long as S ^ [111 |l2j . 

(d) By construction, the bosonic modes of the field 
theory (gj exist at small momenta k s» Hjl |2j| H3 • 
The dynamics at momenta close to tt is captured by the 
superposition of an infinite number of these modes as 
becomes evident from Eq. @ and from Eq. ©. 

The above arguments prove that there can be several, 
rather different looking descriptions for the same physics. 



One may not conclude from the validity of a particular 
description, e.g. in terms of a multitude of modes, that 
another description, e.g. in terms of only a few modes, is 
not valid. 

Bearing these considerations in mind, we turn to a 
continuous unitary transformation (CUT) of the Hamil- 
tonian JQl onto an effective Hamiltonian in terms of 
triplons. To proceed perturbatively we transform into 



H — H' / J — [S 2 iS 2 i+i + AS2iS2i-i] , 



(10) 



where J = J (l + S) and A = (1 - 5)/(l + 5). We 
expand about isolated dimers at A = 0. A CUT is 
used to map the Hamiltonian H to an effective Hamil- 
tonian H e ff which conserves the number of triplets on 
the strong bonds, i.e. [H ,H e g] = where H := H\\= 
(lif. Hence these triplets become the elementary ex- 
citations, triplons, of the effective Hamiltonian. The 
ground state of H c g is the triplon vacuum. An in- 
finitesimal antihermitian generator r\ induces the flow 
dH/dl — [rj(l), H(l)] |2S|. where I is an auxiliary vari- 
able, H(0) = H, H(oo) = H c ff. An appropriate choice 
for rj is 

^■(l) - sgn([Foki - [H ] jd )H itj (l) (11) 

where the matrix elements rjij and H^j are given in an 
eigen-basis of Ho. The choice 111 H retains only triplon 
conserving processes eliminating all parts of H changing 
the number of triplons [l5| . 

In order to determine spectral weights corresponding 
to a given observable R this observable is transformed in 
the same way as the Hamiltonian 



dR/dl = [rj(l),R(l)] 



(12) 



leading for I — > oo to the effective observable i? c ff. Af- 
ter the transformation, i.e. at I = oo, the subspaces with 
different number of triplons are disentangled: the Hamil- 
tonian conserves the number of triplons and does not mix 
subspaces with different numbers of triplons. Given this 
property, it makes sense to define the spectral weight of 
each of these subspaces by 



I n = (0\Rl s P n R ee \0) 



£{0,1,2,3,...} (13) 



where P n projects onto the subspace with n triplons 
U3 E£i uM- This requires only explicit counting 
since the ground state |0) after the transformation is the 
triplon vacuum, i.e. it is the product state of singlets on 
all strong bonds. The technically involved explicit defi- 
nition of "reduced exclusive matrix elements" is avoided 
|29|. Note that the definition of spectral weights like I n 
is useless if the triplon number is not conserved. Only 
the triplon conservation ensures that the total correla- 
tion splits into additive contributions from the various 
subspaces. 
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The CUT cannot be carried out completely. Perturba- 
tively the differential equations are truncated when terms 
of a certain order in A arise 15] . The local spin compo- 
nent Sj is the observable R. The effective observable 
R e g is calculated till order 7 in the one-, two-, and three- 
triplon sector. The four-triplon sector is transformed till 
order 6. 

The total weight I to t = X)o° hi is given by the sum rule 
hot = (R^R) = ((Sj) 2 ) = 1/4 which serves as sensitive 
check for the validity of the results. In the following we 
will discuss the relative weights In,rci = In /hot = 4J n . 
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TABLE I: Relative weight Ii, rc i of the one-triplon channel, 
a) position of the singularity from unbiased approximants; b) 
exponent at the unbiased positions; c) exponents in the biased 
approximants. (blanks: approximants without singularity). 

For isolated dimers (A = 0), the total spectral weight 
lies in the one-triplon channel /i. rc i = 1. As A increases 
ii jrc i decreases and the weight flows into the multi-triplon 
sectors. In Tab. the results for unbiased dlogPade ap- 
proximants indicate a singularity at A ~ 1. From the 
physics of the Hamiltonian JQl we know that the sin- 
gularity is located at A = 1 where the system becomes 
critical. Thus it is advised to investigate approximants 
biased to display the singularity at unity: I\ oc (1 — A) 7 
(last column in Tab. HJ. The exponent is found to be 
7 = 0.325 ± 0.016. which leads us to conjecture that 
it takes exactly the value 7 = 1/3. More generally, we 
presume that any single mode, which vanishes due to 
mixing with a continuum with square-root singularities 
at the band edges, looses its weight like /mode oc V Aw, 
where Aoj is the distance of the mode to the band edge of 
the continuum. The weight in the continuum is assumed 
to be constant. 

We support our claim by considering a generic resol- 
vent I(oj) = l/(w-a-E(o;)) with £(w) = [w+\/w 2 - 4]/2 
for uj < —2. Such a resolvent appears for instance in the 
dynamics of two hard-core particles hopping from site to 
site in one dimension at given total momentum \l2\ . The 
constant a allows to tune a nearest-neighbor interaction 
(a < 0: attraction, a > 0: repulsion) whereas £(w) in- 
corporates the kinetic energy of the relative motion. For 
a < — 1 a bound state emerges from the continuum. It 
is separated from the continuum by the energy Auj — 
-(2 + a+l/a). Its weight J mode is 1/(1 -^£(-2- Aw)) 
implying I moa _ c oc V Aw for Auo — > 0. 

For dimerized spin chains, the single triplon mode is 
separated from the continuum by an energy Aw of the 
order of the energy gap A 01 ■ Using I mo do oc V Alu 



for each total momentum and integrating then over all 
momenta to obtain the local weights we find that the 
single mode looses its weight as vA oc which agrees 
excellently with the extrapolations. 
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TABLE II: Relative weight /2,rci of the two-triplon channel, 
a) position of the singularity from unbiased approximants for 
d\l2,rei; b) exponent at the unbiased positions; c) exponents 
in the biased approximants; d) i2,rei|A=l integrated from the 
biased (position A = 1 and exponent —2/3) approximants (ap- 
proximants without singularity: blanks; with spurious poles: 
stars) . 

In the two-triplon channel, there is no indication for a 
zero of h at A = 1 . On the contrary, Pade approximants 
indicate significant weight at criticality. The weight in 
the two-triplon channel is the weight transferred from 
the one-triplon channel minus the weight transferred fur- 
ther to channels with three and more triplons. Hence it is 
natural to assume the existence of a singularity with ex- 
ponent 1/3 in li. But if this singularity is not linked to a 
zero, dlogPade approximants cannot detect it. Hence we 
investigate the derivative d\Ii which should be governed 
by a divergence with exponent -2/3. Indeed, Pade (not 
shown) and dlogPade approximants (Tab. [Hj indicate a 
singularity at A w 1. Approximants biased to a singular- 
ity of dxh at A = 1 yield exponents 72 ~ —0.76. The 
corresponding value of /2,1'ci found from integrating d\l2 
is 1.25. Since this value overestimates the sum rule by at 
least 25%, we conclude that the exponent 72 ~ —0.76 is 
too large in absolute value. Thus, we bias the approxi- 
mants to the expected behavior d\l2 oc (1 — A)~ 2 / 3 . In 
the last column of Tab.lTTlthe ensuing values for /2,rci are 
given. Quite unexpectedly the results conclusively point 
to a spectral weight very close to unity! Since diagonal 
approximants usually yield the most reliable results we 
retain the value /2.1'ci ~ 0.998, keeping a possible error of 
a few percent in mind. 

The sum rule corroborates the above result strongly. 
Pade and unbiased dlogPade approximands consistently 
show that /3 ilo i is not larger than about 3 • 10~ 4 which 
agrees perfectly with the value close to unity for -?2,rci- 
The inclusion of a singularity does not enhance /3, rc i- 
The biased approximants for d\Iz (position A = 1 and 
exponent —2/3) yield also only a contribution of about 
2 • 10~ 4 ([2,2]). Finally, the Pade and dlogPade approxi- 
mants for the four-triplon contribution Ii iTC i consistently 
indicate values well below 10 -4 . The biased approximant 
for d\h (position A = 1 and exponent —2/3) even yield 
values below 10~ 5 . Therefore, we conclude that the con- 
tributions of channels with four and more triplons can be 
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safely neglected. 

Fig. 1 shows the final results for Ii. rc i, ^2,rel an d ^3,rei; 
the tiny four- and more triplon contributions are ne- 
glected. The sum rule is excellently fulfilled to within 
ps 0.003 for all values of A supporting the above analysis. 



embodies an important message for potential theories for 
high-temperature superconductors. 

We thank M. Gruninger and E. Miiller-Hartmann for 
stimulating and encouraging discussions and the DFG for 
financial support in SP 1073 and in SFB 608. 




FIG. 1: Relative spectral weights J njle i(A) of the dynamic 
structure factor in the dimerized chain. Depicted are biased 
approximants with singularity at A = 1 and exponent 1/3 for 
■fi.rei ([4,2], black solid line); exponent —2/3 in the derivative 
of /2,rci ([3,2], dashed line) and similarly for 1000-/3, re i ([2,2], 
dashed-dotted line). The grey line is the sum X^n=i ^n.rel- 



In conclusion, we have found that even the dynamic 
structure factor of the critical uniform spin chain can 
be described to about 99% by two elementary triplets 
(triplons). This shows that scenario (ii) is correct for 
spin chains. A description in terms of a few triplons is 
possible. This lays the foundation for a new route to 
^spectral densities in many models (first results in |l3l 

Since the spin chain is the archetype of a gapless criti- 
cal model described by spinons we conclude that neither 
the occurrence of fractional excitations nor the vanish- 
ing of the gap precludes the applicability of an approach 
in terms of integer triplons. Even a larger part of the 
spectral weight is covered by two-triplon states than is 
covered by two-spinon states (72.89%) j2||, which calls 
for further investigations of the relation between spinon 
and transformed triplon states. Since these results hold 
in spin chains, the home field of spinons, we conclude 
that a large class of low-dimensional quantum antifer- 
romagnets is accessible by calculations based on integer 
excitations. If this is true for undoped antiferromagnets 
there is no necessity either to resort to fractional exci- 
tations for doped antiferromagnets. Hence our finding 
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